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Warm dark matter is consistent with the observations of the large- 
scale structure, and it can also explain the cored density profiles on smaller 
scales. However, it has been argued that warm dark matter could delay 
the star formation. This does not happen if warm dark matter is made 
up of keV sterile neutrinos, which can decay into X-ray photons and active 
neutrinos. The X-ray photons have a catalytic effect on the formation of 
molecular hydrogen, the essential cooling ingredient in the primordial gas. 
In all the cases we have examined, the overall effect of sterile dark matter 
is to facilitate the cooling of the gas and to reduce the minimal mass of 
the halo prone to collapse. We find that the X-rays from the decay of keV 
sterile neutrinos facilitate the collapse of the gas clouds and the subsequent 
star formation at high redshift. 
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1. Introduction 

Both cold and warm dark matter models agree with the observed struc- 
ture on the large scales. However, there are several inconsistencies between 
the predictions of the cold dark matter (CDM) model and the observa- 
tions [1]. The low cutoff in dark matter contents of dwarf spheroids, the 
smoothness of our dark matter halo, and the old globular clusters (observed 
in Fornax) resisting the infall into the center by dynamical friction [Ij, all 
can be explained by warm dark matter (WDM) because it suppresses the 
structure on scales that are smaller than the free-streaming length. 

While the suppression of the small-scale structure is desirable, it has 
been argued that "generic" WDM (for example, gravitino) can slow down 
structure formation and delay reionization of the universe, which can lead, 
in turn, to an inconsistency with the reionization redshift obtained by the 
WMAP [2j . This problem can be alleviated in the case of the WDM in the 
form of sterile neutrinos with mass of several keV and a small mixing angle 
with the ordinary neutrino O SI [5] because such sterile neutrinos can decay 
and produce photons that catalyze the formation molecular hydrogen and 
speed up the star formation. 

In the absence of metals, gas cooling is mainly due to the collisional ex- 
citation of H2, its subsequent spontaneous de-excitation, and photon emis- 
sion. In the primordial gas clouds, hydrogen molecules can be formed only 
in reactions involving e~ or H"*" as a catalyst. Thus, an X-ray radiation 
can increase the production of the H2 by enhancing the ionization fraction, 
which subsequently leads to speed up of the gas cooling and star formation. 
Although sterile neutrinos are stable on cosmological time scales, they nev- 
ertheless decay. The decay channel important for us is that of decay into 
one active neutrino and one photon, i.e., Vs ^al-, where the photon energy 
is half of the sterile neutrino mass, Eq ms(? jl- These decays produce 
an X-ray background radiation that increases the production of molecular 
hydrogen and can induce a rapid and prompt star formation at high redshift. 

Sterile dark matter has a firm motivation from particle physics [6] [7]. 
The discovery of the neutrino masses implies the existence of right-handed 
gauge-singlet fields, all or some of which can be lighter than the electroweak 
scale. These sterile neutrinos can be produced in the early universe by 
different mechanisms, for example, from neutrino oscillations [6] or from 
the Higgs decays [7], or from the couplings to a low-scale inflaton [8]. The 
same particles, produced in a supernova, could account for the supernova 
asymmetries and the pulsar kicks [9], and can play a role in the formation 
of super-massive black holes in the early universe [lOj. 

We will examine the thermal evolution of the gas clouds, taking into 
account both effects of the sterile neutrino decays, namely, the ionization 
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and heating of the gas. We follow the evolution of the baryonic top-hat 
overdensity, the gas temperature and the H2 and e~ fraction. In order to 
perform the calculation we have incorporated to our previous code |4j, the 
effects of sterile neutrino decays within collapsing halos and absorption of 
the X-ray background from sterile neutrinos by He atoms in the intergalactic 
medium. Our goal is to juxtapose the evolution of the gas temperature in 
the primordial clouds in the CDM model and the WDM model with keV 
sterile neutrinos and estimate of the minimal mass able to collapse at a 
given redshift. 



2. Description of the code 

The top-hat overdensity evolution in a single-zone approximation [IH |1] 
is described by the following equation [12] 



^=^ ("-^"^"^3 -1, (1) 

2(1 — cos a) 

where the parameter a is related to the redshift z and the redshift of viri- 
alization z^ir through 

1 + Zyir _ ^ g - sin g y 
l + z ~\ 27, ) ■ 

According to these equations, the virial value of overdensity {5 ~ ISvr^) is 
reached at the redshift 2:3,r/2 = 1-06555 (1 + z^r) — 1- 

Further evolution of 5 depends on the type of the matter in the over- 
dense region. If it is the dark matter, then after virialization its density 
remains constant forever. The situation is different in the case of baryons. 
If cooling is efficient enough, then the density gradually increases. Other- 
wise, the density remains constant, and there is no star formation in the 
halo. Following [11], we assume that the density of the halo stays constant 
after redshift of -237r/2, i-e., it is equal to the value p = ISvr^OoPo (1 + -^mr) ■, 
where po is the present critical density of the universe. This assumption is 
sufficient for our purposes. 

The evolution of the gas temperature is governed by the equation [3] 



dT , . T dn„ T du T d'j 

— = (7 - 1) ^ + 7 

dz Up dz fi dz (7 — 1) dz 

(7-l)A 

^' 



UpkHo {l + z)JnA + ^o{l + z 
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where T, 7, Up, /i, k and A are the temperature of the gas, the adiabatic in- 
dex, the number density of non-dark matter particles, the molecular weight, 
the Boltzmann constant and the cooling/heating function, respectively. 

If we simply integrate equation ^ to the redshift z^ir then, for primor- 
dial clouds with masses greater than a critical value, the gas temperature 
will be much lower than the virial temperature T^r fT3] 

where fimr and M are the mean molecular weight during virialization and 
the halo mass (combined mass of dark and baryonic constituents). There- 
fore, we must take into account shocks and increase the gas temperature to 
the virial value. We assume that the evolution of the gas temperature is 
linear between the redshifts a^id z^r- 

To calculate the cooling/heating function A in equation ([3]) we follow 
the number density evolution of different baryonic components and take 
into account all relevant chemical and thermal processes, e.g, heating and 
cooling due to the sterile neutrino decays. The detailed list of all included 
processes and their coefficients can be found in [3]. 

Let us define the number fraction of component i SiS Xi — ft i/n, where 
Hi and n = uh + + 2nH2 a-^^ the number density of the ith component 
and hydrogen species, respectively. The time evolution of number fraction 
can be described by the kinetic equation: 

~Tr ~ kimiXlXm + kjiXj, (5) 

I m j 

where the first component on the right-hand side describes the chemical re- 
actions and the other one accounts for photoionization/photodissociation 
processes. Coefficients /c/mi and kji are reaction rates and photoioniza- 
tion/photodissociation rates multiplied by the numbers equal to 0, ±1 or 
±2 depending on the reaction. In our calculations, we have considered 
the following five species: H, H+, H2, e~, and He with the mass fraction 
Y = 0.244 [14J. We used simplified molecular hydrogen chemistry similar 
to the approach presented in [TT] . 

In order to take into account the effects of sterile neutrino decays inside 
the collapsing halo, we have solved the radiative transfer equation for the 
spherically symmetric clouds with uniform density (see Section [3] and also 
[151 116j). In addition, we have included absorption of the X-rays from the 
sterile neutrino decays by both H and He. 

The photons from the decay of the sterile neutrinos are mainly absorbed 
by neutral helium and hydrogen atoms leading to their ionization. The 



main printed on February 2, 2008 



5 



ionization rate due to these photons is enhanced almost 100 times due to 
additional ionization by the secondary electrons, which deposit almost 1/3 of 
their energy into ionization. The energy of the absorbed photons partially 
goes into ionization and partially into heating and excitations. We have 
adopted the approximation [3], in which the ionization rate (in units 
s~^) and heating (in units erg s~^ cm^^) due to the photons from the decay 
of the sterile neutrinos are respectively equal to 



where h, an {i^), hiy^ = 13.6 eV, (u), hv^^ = 24.6 eV and X are the 
Planck constant, photon frequency, the cross section and energy threshold 
for H and He ionization, and hydrogen mass ratio, respectively. The coef- 
ficients d = 0.3908, Oi = 0.4092, bi = 1.7592, Ch = 0.9971, an = 0.2663 
and = 1.3163 are taken from [T7]. The function Aj„t (z) (in units erg 
s^^ cm^^) is the energy absorption rate of the photons from the sterile 
neutrino decays inside the collapsing cloud and is given by equation (p^ . 
Finally, ly (z) (in units of erg cm~^ s~^ sr~^ Hz~^) is the specific intensity 
of the X-ray background from the sterile neutrino decays, which takes into 
account absorption by H and He in the intergalactic medium and which can 
be calculated in a similar way to the specific intensity given in [4J. 



Since the cross section for absorption of the X-ray photons by He is much 
larger than the same cross section for H, a large amount of energy can be 
accumulated in free electrons due to He ionization. Therefore, changes in 
the He ionization can strongly affect the absorbed energy, and we must take 
into account their effects to the heating/cooling function and the ionization 
by the secondary electrons, as in equations ^ and ([7]). 





(7) 
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To estimate the minimal mass of the primordial halo able to collapse at 
a given redshift we use the following criterion of successful collapse [llj: 

T{7]Zyir) < r]T{Zyir) , (8) 

where we take r] = 0.75. It means that the cloud is considered to collapse 
if its temperature drops substantially within a Hubble time, which roughly 
corresponds to the redshift dropping by a factor 2^/^. 



3. Radiative transfer inside the collapsing cloud 

In order to derive the energy absorption rate of the photons from the 
sterile neutrino decays inside the collapsing gas clouds Aint{z), let us as- 
sume that these halos have spherically symmetric shape and their densities 
are uniform. In that case, it will be convenient to consider the radiative 
transport equation in a system of coordinates (r, p) defined by the transfor- 
mation formula: (r, /i) {r,p = r\J\ — for — 1 < /i < 1 [181 115|. where 
r is a radial coordinate, /i = cos0, and Q is the angle between the outward 
normal and the photon direction. 

For a given radius r, the "impact" parameter p can vary between and 
r. Because the parameter p cannot distinguish between fj. > and n < 0, 
the radiation intensity Ii, has to be separated into outward and inward 
/~ directed intensity, respectively. 

Now, we can cast the time-independent, non-relativistic equation for 
radiation transport in spherical geometry into the form 



K{ru,p) (9) 

ju {Tu,P) - Sy , (10) 



where 



jv{r,p) = ^(/+(r,p)+/-(r,p)), (11) 

K{r,p) = l(/+(r,p)-/-(r,p)), (12) 

and dTy = —Xud{rfi), is the optical depth at the radius r. The term 5,^ = 
Vv/Xv is a source function, where r]u (in units of erg s~^ cm~'^ sr^-*^ Hz^^) 
and Xu (in units of cm~^) are the total emissivity and opacity at frequency 
v, respectively. 

We consider sterile neutrino decays inside the spherically symmetric 
cloud, thus, the photon flux from this process will be peaked around the 
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frequency corresponding to the energy of Eq = m<jC^/2. Since the density of 
the cloud is much higher than its surroundings, there wih be only few pho- 
tons at this energy impinging upon the outer boundary, which is set by the 
radius R of the cloud. The external radiation field also consists of the X-ray 
photons emitted due to the sterile neutrino decays at large distances from 
the cloud, however, as they have been emitted at earlier times, they will be 
redshifted to lower energies. Thus, they will not give the contribution to 
the photon flux at the energy of Eq. 

Assuming, that there is no external radiation field at the energy of Eq, 
we can write the boundary conditions as follows [15] 



Kip,p) = 0<p<R, (13) 

jAR,P) = K{R,p) 0<p<R. (14) 

Since we are interested only in the photon flux at the energy of Eq, 
from now on, we drop the z/ dependence for clarity. It means that we have 
to multiply all of the quantities by the Dirac delta S {Eq/K) or understand 
them as they have been already integrated over frequency. We will use the 
latter interpretation. According to the emissivity and opacity definition, we 
have 

_ dn^ Eq _ VtdmC^Q , . 

~ dt iTT- SttTs ' ^ ' 

X = nuCFH {Eq) + nueCTHe {Eq) , (16) 

where g is the total dark and baryonic matter mean density of the collaps- 
ing cloud, is the number density of emitted photons due to the sterile 
neutrino decays. We assume that all of the dark matter in the collapsing 
halo consists of the sterile neutrinos and that the number density of sterile 
neutrinos do not change with time. The latter assumption can be justified 
by the fact that the inverse width of the sterile neutrino radiative decay, 
Ts, is much longer than the age of the universe. In addition, we neglect the 
stimulated emission. 

Equations ([9]) and (jlOp can be rewritten as 



{t,p) 
d'^h {t,p) 



j{r,p)-S, (17) 

-h{T,p), (18) 

which with the boundary conditions ()13p and p4p have the following solu- 
tion 

h{T,p) = -Se^^PhmhT, (19) 
j{T,p) = -Se^(P)coshr + 5, (20) 
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where T (p) = —xRfJ- is the optical depth at the cloud boundary derived 
for the given impact parameter p. The luminosity of the collapsing cloud is 
given by 

L = 16-k'^R^ [ h{R,n)fidn, (21) 
Jo 

whereas if we neglect absorption it would be equal to 

L, = fTT^RS. (22) 

The fraction of the energy absorbed by the cloud is given by 1 — L/Lg, thus 
the energy absorption rate (in units of erg cm~^ s~^) is equal to 

A,„, = (l _ £) (23) 

Using equations ([IS]), ([H]), ([IS]), ([21]) - ([23]), the definition of the source 
function 5, and doing some algebra we get 

A.n*(.) = ^|g^/(a)n(z), (24) 



where 



/(a) = 1-?A 



a (z) 



4 



(25) 



[UH {z) an (Eo) + nne (z) aue {Eq)\ . (26) 

[z) J 



The term mn and M denotes hydrogen mass and the total mass of the 
cloud, respectively. 



4. Results 

We have performed a detailed analysis of the cooling and collapse of 
the primordial gas clouds in the model with warm dark matter, taking into 
account both the increase in the fraction of molecular hydrogen and the 
heating due to the sterile neutrino decays. To illustrate these effects, we 
have performed the analysis for some benchmark cases which arise in real- 
istic scenarios [3]. The effect on the largest gas clouds is negligible, whereas 
smaller clouds will be affected: for the largest clouds, the additional molec- 
ular hydrogen makes no difference, but for the smaller ones, the increase in 
the X-ray background makes the collapse possible in cases where it could 
not occur in the absence of sterile neutrino decays. 
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Fig. 1. Top left, top right and bottom left: Evolution of ionization fraction, H2 
fraction and temperature with redshift for different models. In each case, the 
mass of the primordial cloud is equal to M = lO^M© and virialization redshift 
to Zvir = 12. Bottom right: Dependence of the minimal mass of primordial halo 
able to collapse on its virialization redshift. Models we have used in calculation 
are following: = 25 keV and sin^ 6* = 3 x IQ-^^ (WDMl), = 15 keV and 
s\r?e = 3 X 10-12 (WDM2), = 3.3 keV and sin^ 61 = 3 x 10"^ (WDM3), and 
CDM. 



Our results presented in Fig. [T] show that the overall effect of sterile 
neutrino decays is to enhance the H2 fraction and to speed up the cooling 
of the gas in the primordial halos. The minimal mass of the cloud able to 
collapse is reduced in all WDM models we have examined. We note that 
a more detailed treatment of the sterile neutrino free-streaming may affect 
our results: we did not take into account the filamentary star formation |19j : 
some other effects may also be important j20| . 

In summary, the X-ray photons from sterile neutrino decays could play 
an important role in the formation of the first stars because they increase 
the fraction of molecular hydrogen. 



10 



MAIN PRINTED ON FEBRUARY 2, 2008 



REFERENCES 

[1] X. Hernandez, G. Gilmore, MNRAS 297, 517 (1998); J. Sommer-Larsen, A. 
D. Dolgov, Astrophys. J. 551, 608 (2001); F. Governato et ai, Astrophys. J. 
607, 688 (2004); M. Fellhauer et ai, Astrophys. J. 651, 167 (2006); B. AUgood 
et ai, MNRAS 367, 1781 (2006); T. Goerdt et al., ibid., 368, 1073 (2006); 
G. Gilmore et al., Astrophys. J., 663, 948 (2007); R. F. G. Wyse, G. Gilmore, 
arXiv:0708.1492 [astro-ph]. 

[2] N. Yoshida et al., Astrophys. J. Lett. 591, 1 (2003). 

[3] P. L. Bierniann, A. Kusenko, Phys. Rev. Lett. 96, 091301 (2006). 

[4] J. Stasielak, P. L. Biermann, A. Kusenko, Astrophys. J. 654, 290 (2007). 

[5] J. Stasielak, P. L. Biermann, A. Kusenko, arXiv:astro-ph/0701585, 

[6] S. Dodelson, L. M. Widrow, Phys. Rev. Lett. 72, 17 (1994); K. Abaza- 
jian, G. M. Fuller, M. Patel, Phys. Rev. D 64, 023501 (2001); A. D. Dolgov, 
S. H. Hansen, Astropart. Phys. 16, 339 (2002); X. d. Shi, G. M. Fuller, Phys. 
Rev. Lett. 82, 2832 (1999); T. Asaka, et al., Phys. Lett. B 631, 151 (2005); 
Phys. Lett. 5 620, 17 (2005); Phys. Lett. 5 638, 401 (2006); M. Shaposhnikov, 
Nucl. Phys. B 7 63, 49 (2007); A. Palazzo, D. Cumberbatch, A. Slosar, J. Silk, 
larXiv:0707.T495 [astro-ph]. 

[7] A. Kusenko, Phys. Rev. Lett. 97, 241301 (2006). 

[8] M. Shaposhnikov, I. Tkachev, Phys. Lett. B 639, 414 (2006). 

[9] A. Kusenko, G. Segre, Phys. Lett. B 396, 197 (1997); A. Kusenko, G. Segre, 
Phys. Rev. D59, 061302 (1999); M. Barkovich, J. C. D'Olivo, R. Montemayor, 
Phys. Rev. D 70, 043005 (2004); G. M. Fuller, A. Kusenko, I. Mocioiu, S. Pas- 
coh, Phys. Rev. D 68, 103002 (2003); A. Kusenko, Lnt. J. Mod. Phys. D 13, 
2065 (2004); C. L. Fryer, A. Kusenko, Astrophys. J. Suppl. 163, 335 (2006); 
[10] P. L. Biermann, F. Munyaneza, 'arXiv: astro-ph/0702 1 73, 

[11] M. Tegmark, J. Silk, M. J. Rees, A. Blanchard, T. Abel, F. PaUa, Astrophys. 
.J. 474, 1 (1997). 

[12] T. Padmanabhan, Structure Formation in the Universe, Cambridge University 
Press, Cambridge 1993. 

[13] R. Barkana, A. Loeb, Phys. Rep. 349, 125 (2001). 

[14] Y. I. Izotov, T. X. Thuan, Astrophys. J. 500, 188 (1998). 

[15] J. Stasielak, S. Stachniewicz, M. Kutschera, Acta Phys. Pol. B 36, 3897 
(2005). 

[16] D. Mihalas, B. W. Mihalas, Foundations of Radiation Hydrodynamics, Oxford 
University Press, New York 1984. 

[17] J. M. ShuU, M. E. van Steenberg, Astrophys. J. 298, 268 (1985). 

[18] D. G. Hummer, G. B. Rybicki, MNRAS 152, 1 (1971). 

[19] L. Gao, T. Theuns, Science 317, 1527 (2007). 

[20] E. Ripamonti, M. MapeUi, A. Ferrara, MNRAS 375, 1399 (2007). 



